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Abstract 



We determine by numerical simulations on a lattice the gauge-invariant two-point cor- 
relation function of the gauge field strengths in the QCD vacuum, down to a distance of 
0.1 fm. 
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1. Introduction 



The gauge-invariant two-point correlators of the field strengths in the QCD vacuum are 
defined as 

2W(x) = (0\Tr{G w (x)S(x,0)G 1/a (0)S^x,0)}\0) , (1.1) 

where G w = gT a G® p and T a are the generators of the colour gauge group in the funda- 
mental representation. Moreover, in Eq. (1.1), 

S(x,0) = Pexp [i J dtx^A^xt^j , (1.2) 

with = gT a A a y is the Schwinger phase operator needed to parallel-transport the tensor 
G UCT (0) to the point x. 

These field-strength correlators play an important role in hadron physics. In the 
spectrum of heavy QQ bound states, they govern the effect of the gluon condensate 
on the level splittings |I], 0, They are the basic quantities in models of stochastic 
confinement of colour ||, [5], Q and in the description of high-energy hadron scattering 

A numerical determination of the correlators on lattice (with gauge group SU(3)) 



already exists, in the range of physical distances between 0.4 and 1 fm [11]. In that range 
D mvcr falls off exponentially 

?>np,»<r( x ) ~exp(-|x|/A) , (1.3) 



with a correlation length A ~ 0.22 fm [11 



What makes the determination of the correlators possible on the lattice, with a reason- 
able computing power, is the idea JT2], [OJ of removing the effects of short-range fluctua- 
tions on large distance correlators by a local cooling procedure. Freezing the links of QCD 
configurations one after the other, damps very rapidly the modes of short wavelength, 
but requires a number n of cooling steps proportional to the square of the distance d in 
lattice units to affect modes of wavelength d: 

n~kd 2 . (1.4) 
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Cooling is a kind of diffusion process. If d is sufficiently large, there will be a range of values 
of n in which lattice artefacts due to short-range fluctuations have been removed, without 
touching the physics at distance d; by lattice artefacts we mean statistical fluctuations 
and renormalization effects from lattice to continuum. This removal will show up as a 
plateau in the dependence of the correlators on n. This was the technique successfully 
used in Ref. [ll]. There, the range of distances explored was from from 3-4 up to 7-8 
lattice spacings at (3 ~ 6., which means approximately from 0.4 up to 1 fm in physical 
distance. The lattice size was 16 4 . 

We have now new results on a 32 4 lattice, at (3 between 6.6 and 7.2: at these values 
of j3 the lattice size is still bigger than 1 fm, and therefore safe from infrared artefacts, 
but d = 3,4 lattice spacings now correspond to physical distances of about 0.1 fm. Since 
what matters to our cooling procedure is the distance in lattice units, we obtain in this 
way a determination of the correlators at distances down to 0.1 fm. 



2. Computations and results 

The most general form of the correlator compatible with the invariances of the system is 
1,1,1 

T> W;U(T (x) = (g^9 P a - 9pa9 P u) \D{x 2 ) + V x {x 2 ) ^ 

~\~\X pX u (j p a X pX a (j p V -\- X pX a Cj p V X pXyCj^uj — - . (,2.1 ) 

T> and T>\ are invariant functions of x 2 . We work in the Euclidean region. 

We can define a T>\\(x 2 ) and a T>±(x 2 ) as follows. We go to a reference frame in which 
x M is parallel to one of the coordinate axes, say \i = 0. Then 

1 3 cfD 
P,l = -Y / V 0ifll (x)=V + V 1 +x 2 -^ , 

1 3 

V ± = - J2 V mj {x) =V + V X . (2.2) 

6 %<j=l 

On the lattice we can define a lattice operator , which is proportional to T)p P ^ a in 

the continuum limit, i.e., when the lattice spacing a — >• 0. Since the lattice analogue of the 
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field strength is the plaquette n w (n) (the parallel transport along an elementary square 
of the lattice, lying on the /zp-plane), T>^ a will be defined as 



Vi P ,v*(da) = ^3?{(Tr[n w (n + da)S(n + da,n)Q ua (n)S\n + da,n)))} , (2.3) 
where 3? stands for real part and the lattice operator VL V(T {n) is given by 

n ua (n) = Ul(n) - n w (n) - ~Tr[nt» - rW(n)] . (2.4) 

da is a line parallel to one of the coordinate axes, with integer length d in units of the 
lattice spacing a, joining the two sites to be correlated; S is the parallel transport along 
this line. The inclusion of the operator \£l va (in place of simply putting IT^) on the 
right-hand side of Eq. (2.3) ensures that the disconnected part and the singlet part of the 
correlator are left out. In particular, the subtraction of the singlet part (which is anyway 
a small contribution of order 0(a 12 )) means that we are indeed taking the correlation of 
two operators with the quantum numbers of a colour octet. 

In the naive continuum limit (a — ► 0) we have that 

K P ,uAda) aZo a 4 V^ ua (da) + 0(a e ) . (2.5) 

Making use of the definition (2.1) we can also write, in the same limit, 

Vft(da) ~ a 4 Vn(d 2 a 2 ) + C(a 6 ) , 

11 a— »0 

V L L {da) a ~ o a 4 V ± (d 2 a 2 ) + C(a 6 ) . (2.6) 

Equations (2.5) and (2.6) come from a formal expansion of the operator, and are expected 
to be modified, when the expectation value is computed, by lattice artefacts, i.e., by effects 
due to the ultraviolet cutoff. These effects can be estimated in perturbation theory and 



subtracted |14| . Instead we remove them by freezing the quantum fluctuations at the 
scale of the lattice spacing. After cooling, and T>\ are expected to obey Eq. (2.6). 

We shall omit the details of our cooling procedure, which have been described in 



several papers (see Refs. |I2", [I3|| and references therein). We only remark that it is 



a local procedure, which affects correlations at distances that grow with the number of 
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cooling steps as in a diffusion process. We then expect that, if the distance at which 
we observe the correlation is sufficiently large, lattice artefacts are frozen by cooling long 
before the correlation is affected: this produces a plateau in the correlation versus cooling 
step. Our data are the values of the correlation at the plateau; the error is the typical 
statistical error at the plateau, plus a systematic error which is estimated as the difference 
between neighbouring points at the plateau (whenever the plateau is not long enough, 
looking more like an extremum). Typically the global error is three times larger than the 
statistical error. The typical behaviour of T)» and of T>\ along cooling is shown in Fig. 1. 

We have measured the correlations on a 32 4 lattice at distances ranging from 3 to 14 
lattice spacings and at (3 = 6.6, 6.8, 7.0, 7.2. From renormalization group arguments, 

a = ±-fW ■ (2.7) 

Ax 

As (3 — > oo, f(/3) is given by 

/8 X 51/121 f a \ 

/(/?) = (33 exp (-3^) [1 + 0(1/(3)] , (2.8) 



for gauge group 577(3) and in the absence of quarks. At sufficiently large (3 one expects 
that 



v L ± fwr 4 = a^Ma^ 2 ^) ' (2 ' 9) 



1 ^ ( d 2 

^ u 

where f((3) is given by Eq. (2.8) and terms of higher order in a are negligible. 

In Figs. 2 and 3 we plot respectively V^f{(3)~ 4: and P^/(/3) -4 versus rf p h y s = (dj Al) 
f(f3). In these figures we have also plotted the values of the correlators obtained in Ref. 
flTTf , corresponding to physical distances <i p h ys > 0.4 fm. We have applied a best fit to all 
of these data with the functions 

V(x 2 

V^x 2 ) = Sexpl-ipi + --expl-P] . (2.10) 




5 



We have obtained the following results: 

A h B „ 

^ ~ 3.3 x 10 8 , -jf ~ 0.7 x 10 8 , 

a ~ 0.69 , b ~ 0.46 , 

^ i , A a ~ , 2.11 

A L 182 ' A L 94 ' v ; 

with xV^d.o.f. — 1.7. The continuum lines in Figs. 2 and 3 have been obtained using the 



parameters of this best fit. With the value of Al determined from the string tension [15 
we obtain 

\ A ~ 0.22 fm , A a ~ 0.43 fm . (2.12) 

The correlation length A^, which enters the non-perturbative exponential terms of T> and 
Pi, as well as the magnitude of the coefficients A and B, are compatible with the values 
obtained in Ref. (Tl|. We have also tried a different fit, in which T>\ is described by a 
purely perturbative-like term, while T> is still the sum of a non-perturbative exponential 
term plus a purely perturbative-like term. In other words, we have fixed B = and 
A a = in the parametrization (2.10) for T> and T>\\ 

V(x 2 ) = 

V 1 (x 2 ) = —. ,2.13) 




We have thus obtained: 



-4- 2.7xl0 8 ' 1 3 



Ai ' — A L 183 

a-0.4 , 6-0.3. (2.14) 

The value of xV-^d.o.f. is again acceptable (about 2); the slope of the non-perturbative 
exponential term is practically unchanged. The behaviour for D\\ and V± obtained using 
the parameters from this best fit is close to the continuum lines reported in Figs. 2 and 
3. 

As a final comment we notice that we have been able to observe terms proportional to 
l/|x| 4 in the correlations because we have worked at larger values of (3, where the distance 
between two points (far enough in lattice units so that the correlation is not modified by 
cooling before lattice artefacts are eliminated) is small compared with 1 fm in physical 
units. A larger lattice (32 4 ) has been necessary to avoid infrared artefacts. 
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FIGURE CAPTIONS 



Fig. 1. A typical behaviour of T)k (crosses; d = 6, (3 = 6.6, lattice 32 4 ) and of V>\ (triangles; 
d — 12, (3 — 6.6, lattice 32 4 ) during cooling. 

Fig. 2. The function V^f(P)~ 4 versus physical distance (in fermi units). Crosses correspond 
to (3 = 6.6, triangles to (3 = 6.8, hexagons to (3 = 7.0, diamonds to (3 = 7.2; crossed- 
circles correspond to the data of Ref. [|ll|]. The line is the curve for Dn obtained 
from the best fit of Eqs. (2.10) and (2.11). 

Fig. 3. The function T>^_f(f3)~ 4 versus physical distance (in fermi units). The symbols are 
the same as in Fig. 2. The line is the curve for V± obtained from the best fit of 
Eqs. (2.10) and (2.11). 
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